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Fig. 25. Plasmodium ( Cancer Pagurus ). 

Fig. 26. Corpuscles of Aster acanthion vulgare, freshly drawn. 

Figs. 27—28. Union of a group of corpuscles of Asteracanthion vulgare. 

Fig. 29. Portion of a plasmodium produced by tbe union of the finely granular 
corpuscles of Echinus sphcera , showing distinct endoplasm containing 
the coarsely granular and the coloured corpuscles, and ectoplasm sending 
out filamentous pseudopodia, which unite with those of free corpuscles. 
Figs. 30—33. Phonergates vorax, from “ Zeitsch. f. Wiss. Zool.” Bd. XXX. 1878. 
Taf. II, figs. 54-57. 

All the figures drawn with Yerick, Oc, 2, Obj. 7. 


IV. 44 On the Analytical Expressions which give the History of 
a Fluid Planet of Small Viscosity, attended by a Single 
Satellite/’ By Gr. H. Darwin, F.R.S. Received March 6, 
1880. 

In a series of papers read from time to time during the past two 
years before the Royal Society, I have investigated the theory of the 
tides raised in a rotating viscous spheroid, or planet, by an attendant 
satellite, and have also considered the secular changes in the rotation 
of the planet, and in the revolution of the satellite. Those investi¬ 
gations were intended to be especially applicable to the case of the 
earth and moon, but the friction of the solar tides was found to be a 
factor of importance, so that in a large part of those papers it became 
necessary to conceive the planet as attended by two satellites. 

The differential equations which gave the secular changes in the 
system were rendered very complex by the introduction of solar 
disturbance, and I was unable to integrate them analytically; the 
equations were accordingly treated by a method of numerical quadra¬ 
tures, in which all the data were taken from the earth, moon, and 
sun. This numerical treatment did not permit an insight into all the 
various effects which might result from frictional tides, and an analy¬ 
tical solution, applicable to any planet and satellite, is desirable. 

In the present paper such an analytical solution is found, and is 
interpreted graphically. But the problem is considered from, a point 
of view which is at once more special and more general than that of 
the previous papers. 

The point of view is more general in that the planet may here be 
conceived to have any density and mass whatever, and to be rotating 
with any angular velocity, provided that the ellipticity of figure is not 
large, and that the satellite may have any mass, and may be revolving 
about its planet, either consentaneously with or adversely to the plane¬ 
tary rotation. On the other hand, the problem here considered is 
more special in that the planet is supposed to be a spheroid of fluid of 
small viscosity; that the obliquity of the planet’s equator, the inclina- 
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tion and the eccentricity of satellite’s orbit to the plane of reference are 
treated as being small, and, lastly, it is supposed that the planet is 
only attended by a single satellite. 

The satellite itself is treated as an attractive particle, and the planet 
is supposed to be homogeneous. 

The notation adopted is made to agree as far as possible with that of 
a previous paper, in which the subject was treated from a similarly 
general point of view, but where it was supposed that the equator and 
orbit were co-planar, and the orbit necessarily circular.* 

The motion of the system is referred to the invariable plane, that is, 
to the plane of maximum moment of momentum. 

The following is the notation adopted:—• 

For the planet:— 

M = mass; a — mean radius ; g — mean pure gravity; C = moment 
of inertia (neglecting ellipticity of figure) ; n = angular velocity of 
rotation; i = obliquity of equator to invariable plane, considered as 
small; cj =■ \q\a. 

For the satellite :— 

m = mass; c = mean distance; Q = mean motion ; e = eccentricity 
of orbit, considered as small; j = inclination of orbit, considered as 
small ; t =s m/c 3 , where m is measured in the astronomical unit. 

For both together:—* 

v = M/m , the ratio of the masses; s = f[(ui//gr) 3 (l + 1 /)]^; h — the 
resultant moment of momentum of the whole system ; E — the whole 
energy, both kinetic and potential, of the system. 

By a proper choice of the units of length, mass, and time, the nota¬ 
tion may be considerably simplified. 

Let the unit of length be such that M+m, when measured in the 
astronomical unit, may be equal to unity. 

Let the unit of time be such that s or +may be 

unity. 

Let the unit of mass be such that C, the, planet’s moment of 
inertia, may be unity. 

Then we have 

Q2c 3 =i¥hm=l ....... (1). 


Now, if we put for g its value M/a 2 , and for v its value M/m , we 
have 



a 2 

M _ m / 5 m’ 


since M+m is unity, 


and since s is unity, m ~ f a 2 , wdien m is estimated in the astronomical 
unit. 

* “ Determination of the Secular Effects of Tidal Friction by a Graphical Me¬ 
thod,” “ Proc. Roy. Soc.,” No. 197, 1879. 
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Again, since C = -f-Mh 2 , and since C is unity, therefore A/==f/a 3 , where 
M is estimated in the mass unit. 

Therefore is unity, when M and m are estimated in 

the mass unit, with the proposed units of length, time, and mass. 

According to the theory of elliptic motion, the moment of momentum 
of the orbital motion of the planet and satellite about their common 

centre of inertia is -tit: —Oc 2/ ^l — e 2 . ISTow it has been shown that the 

factor involving M is unity, and by (1) Qc 2 —Q~%=c%. 

Hence, if we neglect the square of the eccentricity e, the m. of m. of 
orbital motion is numerically equal to Q"a or ch 

Let x=Q ~a=ch 

In this paper x, the moment of momentum of orbital motion, will 
be taken as the independent variable. In interpreting the figures 
given below it will be useful to remember that it is also equal to the 
square root of the mean distance. 

The moment of momentum of the planet’s rotation is equal to Cn ; 
and since 0 is unity, n will be either the m. of m. of the planet’s rotation, 
or the angular velocity of rotation itself. 

With the proposed units r =fm/c 3 =fa< 2 af~ 6 , since m— fa 2 ; and 
a—%Mj a 3 =fm. 

Also T 2 /g (a quantity which occurs below) is equal to cfi/vx 12 . 

Now let t be the time, and let 2 f be the phase-retardation of the 
tide which I have elsewhere called the sidereal semi-diurnal tide of 
speed 2 n 9 which tide is known in the British Association Iteport on 
Tides as the faster of the two K tides. 

Then if the planet be a fluid of small viscosity, the following are 
the differential equations which give the secular changes in the 
elements of the system: 


dn _ 

dt 

dx _ 

dt 

dt 

dt 

1 de _ 

e dt 


. to - 

i£sin4/(l-£).(3). 

Q \ n J 

i^sin4/(i±i)(l- 2n ) .... (4). 

— p JL. sin 4/ -C-~t 'll ........ ( 5 ). 

S « 

am 4>f.U ..... ( 6 ). 

g * x\ n ) 


The first three of these equations are in effect established in my 
paper on the “Precessionof a Viscous Spheroid,§ 17, p. 497, eq. (80). 

* “ Phil. Trans.,” Part II, 1879. 
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The suffix m 3 to the symbols i and N there indicates that the equa¬ 
tions (80) only refer to the action of the moon, and as here we only 
have a single satellite, they are the complete equations. N is equal to 
n/n 0 , so that n 0 disappears from the first and second of (80) ; also 
lUL—l/snQQff, and thus n 0 disappears from the third equation. P = cos?‘, 
Q =sin i, and, since we are treating i the obliquity as small, P=l, 
Q—i ; also \—Q/n\ the e of that paper is identical with the / of the 
present one; lastly £ is equal to Q ( /Q~3, and since with our present 
units s— 1, therefore jmd^/dt=dQ~^ln 0 dt~dxln 0 dt. 

With regard to the transformation of the first of (80) into (4) of 
the present paper, I remark that treating i as small ^PQ— -|\Q 
— 4^(1 — 20/%), and introducing this transformation into the first of 
(80), equation (4) is obtained, except that i occurs in place of (i -f j ). 
Now in the paper on the “ Precession of a Viscous Spheroid ” the in¬ 
clination of the orbit of the satellite to the plane of reference was 
treated as zero, and hence j was zero ; but I have proved in a paper 
“ On the Secular Changes in the Elements of the Orbit of a Satellite 
revolving about a tidally distorted Planet” (read before the Boyal 
Society on December 18th, 1879, but as yet unpublished) that when 
we take into, account the inclination of the orbit of the satellite, the P 
and Q on the right-hand sides of eq. (80) of “ Precession ” must be 
taken as the cosine and sine of i-\-j instead of i. Equations (5) and 
(6) are proved in § 10, Part II, and § 25, Part V of the unpublished 
paper, and the reader is requested to take them as established. 

The integrals of this system of equations will give the secular 
changes in the motion of the system under the influence of the 
frictional tides. The object of the present paper is to find an analy¬ 
tical expression for the solution, and to interpret that solution geome¬ 
trically. 

From equations (2) and (4) we have 


But from (3) and (5) xdjjdt+j dxjdt is equal to the same expres¬ 
sion; hence 


. dn di dj . .dx 

% -j- n —— x — pi —. 

dt dt dt . dt 


The integral of this equation is in—jx. 


or 


i _ x 

j n 


(SO- 


Equation (7) may also be obtained by the principle of conservation 
of moment of momentum. The motion is referred to the invariable 
plane of the system, and however the planet and satellite may interact 
on one another, the resultant m. of m. must remain constant in 
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direction and magnitude. Hence if we draw a parallelogram of 
which, the diagonal is h (the resultant m. of m. of the system), and of 
which the sides are n and as, inclined respectively to the diagonal at 
the angles i and j, we see at once that 

sin?!_ x 

sin j n 

If i and j be treated as small this reduces to (7). 

Again the consideration of this parallelogram shows that 

h a? 3 + 2nx cos (?! +j ), 

which expresses the constancy of moment of momentum. If the 
squares and higher powers of i+j be neglected, this becomes 

}i—n-\~x .(8). 


Equation (8) may also be obtained by observing that dnjdt+dxldt— 0, 
and therefore on integration n + x is constant. It is obvious from the 
principle of m. of m. that the planet’s equator and the plane of the 
satellite’s orbit have a common node on the invariable plane of the 
system. 

If we divide equations (4) and (6) by (8), we have the following 
results:— 


1 dK = 1 2Q 

i dx 2 n\ i) n —Q 


(9). 


lde_l lln—180 

e dx 2x n —Q 


( 10 ). 


But from (7) and (8) 

i+i=i+-=-, 

i x x 


also Q = ar 3 , and n—h—x. 

Hence (9) and (10) may be written 


How 

Therefore 


d_ IQ h 1 x*Qi-x)- 2 \ 

dx ^ 2 x(h — x) x^(h — x) — 1 I 

cl \ = JL ( 

dx °^ 6 2x x s (h—x) — 1 / 

h{x*( h-x)- 2} _ h h x* ' 

2x(h—x){x s (h—x) — l} x(h—x) 2 x^—hx^ + l 

* 

d -i •_1 | 1 ■ li/ x 

^iog 

!W(h-x)-l8 _9_ 7 x*(x~K) 

*2®{* 3 (7t—a)—1} * 2 i6 4 —7i® 3 + l 


. . ( 11 ). 


( 12 ). 


Also 
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cl , 9 7 x 2 (x-h) 

- - log 1 e =~—— .--- -- 

dx x 2 a? 4 — hx*-\- L 


[Mar. 18, 
• (13). 


These two equations are integrable as they stand, except as regards 
the last term in each of them. 

It was shown in a previous paper that the whole energy of the 
system, both kinetic and potential, was equal to \ [n 2 —af 2 ].* 

Then integrating (12) and (13), and writing down (7) and (8) 
again, and the expression for the energy, we have the following equa¬ 
tions, which give the variations of the elements of the system in terms 
of the square root of the satellite’s distance, and independently of the 
time. 


log f 


’ h—x 


+ 


§7l ; 


xhlx 


log 6 — log X® — ~ 


ll — X . 


x 4 —hx $-\-1 

+ const. 


+ const. 


fx 2 (x—li)dx 

jad-^ + 1 


(14). 


n—h — x. 

2E=(h-xy--~-. 


When the integration of these equations is completed, we shall, 
have the means of tracing the history of a fluid planet of small 
viscosity, attended by a single satellite, when the system is started 
with any given moment of momentum and with any mean distance 
and (small) inclination and (small) eccentricity of the satellite’s 
orbit, and (small) obliquity of the planet’s equator. It may be re¬ 
marked that h is to be taken as essentially positive, because the sign 
of h merely depends on the convention which we choose to adopt as 
to positive and negative rotations. 

These equations do not involve the time, but it will be shown later 
how the time may be also found as a function of x. It is not, how¬ 
ever, necessary to find the expression for the time in order to know 
the sequence of events in the history of the system. 

Since the fluid which forms the planet is subject to friction, there¬ 
fore the system is non-conservative of energy, and therefore x must 
change in such a way that JE may diminish. 

If the expression for E be illustrated by a curve in which E is the 
vertical ordinate and x the horizontal abscissa, then any point on this 
“ curve of energy ” may be taken to represent one configuration of 
the system, as far as regards the mean distance of the satellite, 

* c< The Secular Effects of Tidal Friction/’ k\, “ Proc. Roy Soc./’ Mo. 197, 
18/9, eq. (4). 
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Then such a point must always slide down a slope of energy, and we 
shall see which way x must vary for any given configuration. This 
consideration will enable us to determine the sequence of events, 
when we come to consider the expressions for i, e, j, n in terms of x. 

We have now to consider the further steps towards the complete 
solution of the problem. 

The only difficulty remaining is the integration of the two expressions 
in the first and second of (14). From the forms of the expressions 
to be integrated, it is clear that they must be split up into partial 
fractions. The forms which these fractions will assume will of course 
depend on the nature of the roots of the equation — hx ^-\-1 = 0. 

Some of the properties of this biquadratic were discussed in a 
previous paper, but it will now be necessary to consider the subject 
in more detail. 

It will be found by Ferrari’s method that 

ha? + 1 = (» 3 +2a£h- + 1 ( - 2a£?±Z' + M+Z* 1 
l 4 2X.4 / l 4 2V / 

where X 8 --4X—-/i 3 =0. 

By using the property (Xt— li) (Xt + A)=4X, this expression may be 
written in the form 

[{»+i(X*-A)} B + {i(X*-A) v / l+2/a-5}2] x 

[{*-i(x*+A)} 8 + 

Which is of course equivalent to finding all the roots of the bi¬ 
quadratic in terms of h and X. 

Itfow let a curve be drawn of which lv j is the ordinate (negative 


Fig. l. 



The ordinates are drawn to one-third of the scale to which the abscissa? are drawn.. 
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values of 1$ being’admissible) and X tbe abscissa ; it is shown in fig. 1. 
Its equation is & a =X(\ 2 —4). 

It is obvious that OA=OA'=2. 

The maximum and minimum values of W J (viz., B b, B’b ') are given 
by 3X 3 =4 or X=+2/3*. 

Then Bb =BA' = - 2 s /3*+4.2/3* = (4/3*) 3 . 

Since in the cubic, on which the solution of the biquadratic depends, 
W is necessarily positive, it follows that if h be greater than 4/3* the 
cubic has one real positive root greater than OM, and if h be less than 
4/3% it has two real negative roots lying between 0 and OA', and 
one real positive root lying between OA and OM. 

To find OM we observe that since 7i 3 is equal to (4/3^) 2 , and since 
the root of X s •—4X— h 9j ~0 which is equal to —2/3* is repeated twice, 
therefore, if e be the third root (or OM) we must have 

(x+-|) , (x-.)=x»-*-(*) s , 

whence (2/3*) 3 e= (4/3 5 ) 3 , and e or OM=4/3*. 

How OA=2 ; hence, if h be less than 4/3% the cubic has a positive 
root between 2 and 4/3% and if h be greater than 4/3 1 , the cubic has a 
positive root between 4/3* and infinity. 

It will only be necessary to consider the positive root of the cubic. 
How suppose h to be greater than 4/31 

Then it has just been shown that X is greater than 4/3% and hence 
(X being positive) 3\ 3 is greater than 16\, or 4(X 8 —4X) greater than 
X s , or 47A greater than X 3 , or 2 h\~% greater than unity. 

Therefore {^(Xf h) */\ —27A~t} 3 — — {^(Xt + Ji) 27A“£—1 } 3 . 

Thus the biquadratic has two real roots, which we may call a and b. 

Then a=p+i.)[l+v / irn.], 

h=i(\l + 70[l- v^Xri-l]. 

It will now be proved that a is greater and b less than %h. 

How a> or <|7% 

as (Xf + h) [1 -fi a/ 2/A~l—1] > or <37% 

(Xi -j- Jl') , - 

as ■— —| > or < 2A-XI, 

asXI+7i> or < \*^2h —XI, 
as X 3 -f 27zXt+A 3 > or <2fcXi—X 8 , 
as 2x 3 -f 1$ > or < 0. 

Since the left hand side is essentially positive, a is greater than |A. 
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Again b> or < f 7g 

as (Al + 7i)[l— ^ 27?,X“f—1] > or < 37g 
(Xt -}- hi) / -— 

as — —p— v 2 h —Xt > or <2Ji —XI 

Since tlie left-hand side is negative and the right positive, the left is 
less than the right, and therefore b is less than fh. 

If, therefore, h be greater than 4/3^, we may write 

as 4 — hx*+\ — (x— a) (as—b)[(as—«) 3 + /3 3 ], 

where a—f7g —b are positive, and where a is negative. 

We now turn to the other case and suppose h less than 4/3*. All 
the roots of the biquadratic are now imaginary, and we may put 

X^ — JlX* + 1 = [ (x — a) 3 + /3 3 ] [ (x — q) 3 + £ 3 ] . 

If a be taken as — J(Xt — li) 7 then 7 is J(Xt + &). 

Then it only remains to prove that 7 is greater than f h. 

Now 7 > or <f-7g 

as XI > or < 2 h, 
asX 3 > or < 4& 3 =4(X 3 —4X), 
as 16 > or < 3X 3 , 
as 4/3^ > or <X, 

but it has been already shown that in this case, X is less than 4/3^ 
wherefore 7 is greater than f h. 


We may now proceed to the required integrations. 

First case where h is greater than 4/3*. 

Let x 4 —hx s -j-l = (x— a) (as—b)[(as—a ) 3 + j 8 3 ], 

so that the roots are a, b, ajrfit. 

Also let a be the root which is greater than f h, b that which is less, 
and let 

a=a 1 -fj^, b=|7i—b 1? a=%h—cc 1 . 

To find the expression for i we have to integrate — 7 

r b as 4 — hx % +l 

Let f(x) = (x~ a)^($), and let a? 3 //(a?) = Aj(x— a) + 0(a3)/yr(o?). 

Then a ? 3 (x —a) = Af(x) + (x— a) 3 0 (x). 

Hence A —a 3 //' (a). 

If, therefore, /(aj)=® 4 —7&b 3 + 1, A=l/(4a—3fe)=l/4a 1 . 

t 2 
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Thais the partial fractions corresponding to the roots a and b are 

^ 1 __!_ 1 . (15) 

4a 1 aj—a 4^0?—b 


If the pair of fractions corresponding to the roots + be formed 
and added together, we find 

_ i #1(3? ot) -j- n 

2 («i* + / 8 ») [(®-*) 8 + i 8 »] ’ ' ' ’ ^ ^ 

The sum of (15) and (16) is equal to ~ 4 —^TjTi’ an< ^ 


a? 4 —/lar’-j-l 4a 


I log(*cv> a) J-l°g(® w b) log[(*-«*)2 +/3 3 ] 


4(ap + /3 3 ) 


/5 , £U —a 

J-arc tan —_ 


K«j*+P) 5 

Substituting in the first of (14) we have 

(oj co a)8al exp. r.——j——arc tan E 
i=zA—— . _ L4(ai 2 + /3 3 ) _ 

h —a? h /w i 

(# Cob)8bi [(& — a) 3 -}- /5 3 ] 8(a x 9 +/33) 


ft /3 , X — ct 

--arc tan- 

2 + /3 2 ) _A_: 

ho.± 


■ (18). 


where 4 is a constant to be determined by the value of i, which 
corresponds with a particular value of x. 

From the third of (14) we see that by omitting the factor 
xl(Ji—x ) from the above, we obtain the expression for j. 

To find the expression for e we have to integrate .jLJEl ^L. 

F b a 4 -fta 3 -fl 


Now =i(4® 8 -3^) -ifta? 3 , 

and therefore 

The integral remaining on the right hand has been already deter- 
mined in (17). Then substituting in the second of (14), we have 


Bx» 

x*-ka? + 1 )«- 


[.t^V 3f 


• (19)- 


(^C0b) 81j i \_(X — a) 3 -f- / 3 3 ] 8( ai a +/3 = ) 


where B is a constant to be determined by the value of e, correspond¬ 
ing to some particular value of x. 
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From this equation we get the curious relationship 

_B 03 9 

Ai ■(**-&»»+ 1)^ 4 


( 20 ). 


This last result will obviously be equally true even if all the roots 
of a? 4, — hx? + 1 = 0 are imaginary. 

In the present case the complete solution of the problem is com¬ 
prised in the following equations :— 


i=A 


Jl f x—c 

(x cc a)8ai exp. arc tan ~J~ 


A hai 

(X(X) b)8bi [($— Oj) 2 + ^ 3 ]8(a 1 3 +3 B ) 


c _JB a ? 9 
Al(x±-hx s + 

n—h—x. 

2E=(h-xy-\. 

X" 



It is obvious that the system can never degrade in such a way 
that x should pass through one of the roots of the biquadratic 
;A—/i»3 + i == o. Hence the solution is divided into three fields, viz., 
(i) a? =+oo to $=a; here we must write x —a, a?— b for the $ 00 a, 
a?cob in the above solution; (ii) aj=atoas=b; here we must write 
a —a?, a?—b (this is the part which has most interest in application to 
actual planets and satellites) ; (iii) x~~b to $= — go ; here we must 
write a— x, b—x. When x is negative the physical meaning is that 
the revolution of the satellite is adverse to the planet’s rotation. 

By referring to (4) and ( 6 ), we see that i must be a maximum or 
minimum when ?z= 2 Q, and e a maximum or minimum when ^=|fO. 
Hence the corresponding values of x are the roots of the equations 
af*—kx 8 + 2 = 0 , and x*—hx i + - H —0 respectively. 

Since 

it? __ 1 1 11 , 1 — a, Y (x— cc") + / 1 3 

a? 4 — ha? +1 4a x x —a 4 ^ 0 ?—b 2 (a 1 3 + ^ 2 ) [(&—«) 3 + /i 3 ] 
Therefore 

4aj 

-4 --—-— j” — x-i(x —a) H - —a) (-0 — b). 

r 2 (ai a +^) L 1V ' 
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Hence the coefficient of on the right-hand side nrnst be zero, and 


therefore 


_ 1 _ 

4ai 


JL—__f!_=0. 

4bi 2 (»!*+/?) 


A , h lb , 7^ 

And 

ISTow when a?= + oo, arc tan -|?r, and when 03—oo, it is 

equal to —-J-t r. 

Hence when a;=:+ao, ^=J. exp. [ + w7&/}/8(a 1 3 +/} a )], i=—j ; the 
upper sign being taken for + oo and the lower for — oo. 

Then since j tends to become constant when x= + 00 ? and since 
9—7~ o, therefore when a? is very large e tends to vary as x V. 

If x be very small j has a finite value, and i varies as x, and e varies 
as a? 9 . 

j , i, and e all become infinite when x — b, and 1 also becomes infinite 
when x=h. 

This analytical solution is so complex that it is not easy to under¬ 
stand its physical meaning; a geometrical illustration will, however, 
make it intelligible. 

The method adopted for this end is to draw a series of curves, the 
points on which have x as abscissa and i, j , e, n, ]]] as ordinates. The 
figure would hardly be intelligible if all the curves were drawn at once, 
and therefore a separate figure is drawn for i, j , and e ; but in each 
figure the straight line which, represents n is drawn, and the energy 
curve is also introduced in order to determine which way the figure is 
to be read. The zero of energy is of course arbitrary, and therefore 
the origin of the energy curve is in each case shifted along the vertical 
axis, in such a way that the energy curve may clash as little as- 
possible with the others. 

It is not very easy to select a value of h which, shall be suitable for 
drawing these curves within a moderate compass, but after some con¬ 
sideration I chose h— 2'6, and figs. 2, 3, and 4 are drawn to illustrate 
this value of h. If the cubic X 3 —4\— (2*6) 3 =0, be solved by Cardan’s 
method, it will be found that X=2*5741, and using this value in the 
formula for the roots of the biquadratic we have 


$4—2*6 +1 = (a?- 2*539) (aj--826)[(as + *3.82)3 + (-575) 2 ], 

Hence a = 2‘539, b=*826, *= —-382, 0=*575, |7t=r95, and 
4*4=2-356, 414=4-496, ^=2-332, *^ + 0*=5 771. 

Then we have 
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._ , (2’539cvd*)' 553 exp. [’062 arc tan (l’740i»-f ’665) j 

J ~ J ' (®<x> -826)'* 89 (® 2 + 765® + -477)' 18M 


_ B a; 9 . 7 } 

6 At (* 4 —2‘6a; 3 + l)» 

n~ 2*6—a?. 

2i?=(2-6-as) 3 -—h 

ar / 


( 22 ). 


The maximum and minimum values of i are given by the roots of 
tbe equations 4 —2*6a5 3 + 2=0, viz., x = 2‘46 7 and $=1*103. The max¬ 
imum and minimum values of e are given by the roots of the equation 
$ 4 —2*6;r-f xf=0, viz., a?=2*495 and 05=1*0095. The horizontal asymp¬ 
totes for i/A and j/A are at distances from the axis of x equal to 
exp. (*062 x^7t) and exp. (—*062 x^7t), which are equal to 1*102 and 
*908 respectively. 

Fig. 2 shows the curve illustrating the changes of i, the obliquity of 
the equator to the invariable plane. 

The asymptotes are indicated by broken lines; that at A is given 
by x =*826, and is the ordinate of maximum energy; that at B is given 
by a5=2*6, and gives the configuration of the system for which the 
planet has no rotation. The point C is given by $=2*539, and lies on 
the ordinate of minimum energy. Geometrically the curve is divided 
into three parts by the vertical asymptotes, but it is further divided 
physically. 

The curve of energy has four slopes, and since the energy must 
degrade, there are four methods in which the system may change, 
according to the way in which it was started. The arrows marked on 
the curve of obliquity show the direction in which the curve must be 
read. 

Since none of these four methods can ever pass into another, this 
figure really contains four figures, and the following parts of the figure 
are quite independent of one another, viz.: (i) from — oo to O ; (ii) 
from A to O ; (iii) from A to C ; (iv) from + oo to C. The figures 
3 and 4 are similarly in reality four figures combined. For each of 
these parts the constant A. must be chosen with appropriate sign; but 
in order to permit the curves in fig. 2 to be geometrically continuous 
the obliquity is allowed to change sign. 

The actual numerical interpretation of this figure depends on the 
value of A. Thus if for any value of x in any of the four fields the 
obliquity has an assigned value, then the ordinate corresponding to 
that value of x will give a scale of obliquity from which all the other 
ordinates within that field may be estimated. 
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Fig-. 2. 



Diagram for Obliquity of Planet’s Equator.—First case. 

As a special example of this we see that, if the obliquity be zero at 
any point, a consideration of the curve will determine whether zero 
obliquity be dynamically stable or not; for if the arrows on the curve 
of obliquity be approaching the axis of x , zero obliquity is dynamically 
stable, and if receding from the axis of x, dynamically unstable. 

Hence from x— 4- go to B, zero obliquity is dynamically unstable, 
from — oo to 0 and A to 0 dynamically stable, and from A to B, first 
stable, then unstable, and finally stable. 

The infinite value of the obliquity at the point B has a peculiar 
significance, for at B the planet has no rotation, and being thus free 
from what Sir William Thomson calls “gyroscopic domination,” the 
obliquity changes with infinite ease. In fact at B the term equator 
loses its meaning. The infinite value at A has a different meaning. 
The configuration A is one of maximum energy and of dynamical 
equilibrium, but is unstable as regards mean distance and planetary 
rotation ; at this point the system changes infinitely slowly as regards 
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time, and therefore the infinite value of the obliquity does not indicate 
an infinite rate of change of obliquity. In fact if we put n—Q in (1) 
we see that dijidt— —]r(T 3 /g) sin 4/. However, to consider this case 
adequately we should have to take into account the obliquity in the 
equations for dn/dt and dxjdt , because the principal semi-diurnal tide 
vanishes when w=Q. 

Similarly at the minimum of energy the system changes infinitely 
slowly, and thus the obliquity would take an infinite time to vanish. 

We may now state the physical meaning of fig. 2, and this interpre¬ 
tation may be compared with a similar interpretation in the paper on 
u The secular effects of tidal friction,” above referred to. 

A fluid planet of small viscosity is attended by a single satellite, and 
the system is started with an amount of positive moment of momentum 
which is greater than 4/3% with our present units of length, mass and 
time. 

The part of the figure on the negative side of the origin indicates a 
negative revolution of the satellite and a positive rotation of the 
planet, but the m. of m. of planetary rotation is greater (by an amount 
li) than the m. of m. of orbital motion. Then the satellite approaches 
the planet and ultimately falls into it, and the obliquity always 
diminishes slowly. The part from O to A indicates positive rotation 
of both parts of the system, but the satellite is very close to the planet 
•and revolves round the planet quicker than the planet rotates, as in 
the case of the inner satellite of Mars. Here again the satellite 
approaches and ultimately falls in, and the obliquity always diminishes. 

The part from A to C indicates positive rotation of both parts, but 
the satellite revolves slower than the planet rotates. This is the case 
which has most interest for application to the solar system. The 
•satellite recedes from the planet, and the system ceases its changes 
when the satellite and planet revolve slowly as parts of a rigid body— 
that is to say, when the energy is a minimum. The obliquity first 
decreases, then increases to a maximum, and ultimately decreases to 
zero.* 

The part from infinity to C indicates a positive revolution of the 
satellite, and from infinity to B a negative rotation of the planet, but 
from B to C a positive rotation of the planet, which is slower than the 
revolution of the satellite. In either of these cases the satellite 
approaches the planet, but the changes cease when the satellite and 
planet move slowly round as parts of a rigid body—that is to say, when 
the energy is a minimum. If the rotation of the planet be positive, 
the obliquity diminishes, if negative it increases. If the rotation of 
the planet be nil , the term obliquity ceases to have any meaning, since 
there is no longer an equator. 

* According to the present theory, the moon, considered as being attended by the 
.earth as a satellite, has gone through these changes. 
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Pig. 3 illustrates tlie changes of inclination of the satellite’s orbit, 
and may be interpreted in the same way as fig. 2. It appears from the 
part of the figure for which w is negative, that if the revolution of the 
satellite be negative, and the rotation of the planet positive, but the 


Fig. 3. 



Diagram for Inclination of Satellite’s Orbit.—First case. 


m. of in. of planetary rotation greater than that of orbital motion,, 
then, as the satellite approaches the planet, the inclination of the orbit 
increases, or zero inclination is dynamically unstable. In every other 
case the inclination will decrease, or zero inclination is dynamically 
stable. 

This result undergoes an important modification when a second 
satellite is introduced, as will appear in the unpublished paper. 

Pig. 4 shows a similar curve for the eccentricity of the orbit. Tho 
variations of the eccentricity are very much larger than than those of 
the obliquity and inclination, so that it was here necessary to draw the 
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ordinates on a much reduced scale. It was not possible to extend the 
figure far in either direction, because for large values of x, e varies as 
a high power of x (viz., V 1 )* The curve presents a resemblance to that 
of obliquity, for in the field comprised between the two roots of the 
biquadratic (viz., between A and 0) the eccentricity diminishes to a 
minimum, increases to a maximum, and ultimately vanishes at CL 


Fig. 4. 



Diagram for Eccentricity of Satellite’s Orbit.—First case. 


This field represents a positive rotation both of the planet and satel¬ 
lite, but the satellite revolves slower than the planet rotates. This part 
represents the degradation of the system from the configuration of 
maximum energy to that of minimum energy, and the satellite recedes 
from the planet, until the two move round slowly like the parts of a 
rigid body. 

In every other case the eccentricity degrades rapidly, whilst the 
satellite approaches the planet. 

The very rapid rate of variation of the eccentricity, compared with 
that of the obliquity would lead one to expect that the eccentricity of 
the orbit of a satellite should become very large in the course of its 
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evolution, whilst the obliquity should not increase to any very large 
-extent. But it must be remembered that we are here only treating a 
planet of small viscosity, and it will appear, in the unpublished paper 
above referred to, that the rate of increase or diminution of the eccen¬ 
tricity is very much less rapid (per unit increase of x ) if the viscosity 
be not small, whilst the rate of increase or diminution of obliquity 
(per unit increase of x) is slightly increased with increase of viscosity. 
Thus the observed eccentricities of the orbits of satellites and of 
obliquities of their planets cannot be said to agree in amount with the 
theory that the planets were primitively fluids of small viscosity, 
though I believe they do agree with the theory that the planets were 
fluids or quasi-solids of large viscosity. 


We now come to the second case , where li is less than 4 / 3 *. The 
biquadratic having no real roots, we may put 

ad—-f-1 = [ (x — a ) 2 4 / 3 3 ] [ (x —7) 3 4^ 3 ]. 

It has already been shown that a is negative, and 7 greater than f h. 
Let a—%h — a v 7 = 7 1 + 

Then by inspection of the integral in the first case we see that 
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The rest of equations ( 21 ), which express the other elements in 
terms of j and a?, remain the same as before. 

By comparison with the first case, we see that 
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On multiplying both sides of this identity by ad— hx s 41 ? and equating 
the coefficients of ad, we find 
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tlie upper sign being taken for +oo, and the lower for — oo . This 
expression gives the horizontal asymptotes for j and i. 

In order to illustrate this solution, I chose h— 1, and found by tri¬ 
gonometrical solution of the cubic X 8 —4X—1=0, X—21149, and 
thence 

-2-038®+ l-401\ow_ r . Aai 


, a*+ 1*038®+ 714, 


exp. [*081 arc tan (1*500® -f *778) 
+ * 846 arc tan (1*659®—1 * 691) ]. 


B_ _ ^ _ 

a* (® 4 —® 3 +i 

n—l—x. 


2J0=(l-a)*-4 

When x— -f oo, ^7^=1*956= —i/A, and 
when &= oo, y/^l= , 512= — i/A. 



Diagram for Obliquity of Planet’s Equator.—Second case. 

These solutions are illustrated as in the previous case by the three 
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figures 5, 6, 7. There are here only two slopes of energy, and hence 
these figures each of them only contain two separate figures. 

Fig. 5 illustrates the changes of i, the obliquity of the equator to 
the invariable place. 

In this figure there is only one vertical asymptote, viz., that corre¬ 
sponding to oj= 1. For this value of x the planet has no rotation, is 
free from “gyroscopic domination,” and the term equator loses its 
meaning. 

The figure shows that if the rotation of the planet be negative, but 
the m. of m. of planetary rotation less than that of orbital motion, 
then the obliquity increases, whilst the satellite approaches the planet. 

This increase of obliquity only continues so long as the rotation of 
the planet is negative. The rotation becomes positive after a time, and 
the obliquity then diminishes, whilst the satellite falls into the planet. 
In the corresponding part of fig. 2 the satellite did not fall into the 
planet, but the two finally moved slowly round together as the parts of 
a rigid body. 

If the revolution of the satellite be negative, and the rotation of 
the planet positive, but the m. of m. of rotation greater than that of 
revolution, the obliquity always diminishes as the "satellite falls in to 
the planet. 


Fia. 6. 



Diagram for Inclination of Satellite’s Orbit.—Second case. 









275 


1880.] History of Planet and Single Satellite . 

Figs. 2 and 5 only differ in the fact that in the one there is a true 
maximum and a true minimum of obliquity and energy, and in the 
other there is not so. In fact, if we annihilate the part between the 
vertical asymptotes of fig. 2 we get fig. 5. 

Fig. 6 illustrates the changes of inclination of the orbit. It does 
not possess very much interest, since it simply shows that however the 
system be started with positive revolution of the satellite, whether the 
rotation of the planet be positive or not, the inclination of the orbit 
slightly diminishes as the satellite falls in. 

And however the system be started with negative revolution of the 
satellite, and therefore necessarily positive rotation of the planet, the 
inclination of the orbit slightly increases. Fig. 6 again corresponds 
to fig. 3, if in the latter the part lying between the maximum and 
minimum of energy be annihilated. 


Fig, 7. 



Diagram for Eccentricity of Satellite’s Orbit.—Second case. 


Fig. 7 illustrates the changes of eccentricity, and shows that it 
always diminishes rapidly however the system is started, as the satellite 
falls into the planet. This figure again corresponds with fig. 4, if in 
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the latter-i-lie parts between the maximum and mininmm of energy be 
annihilated. 

These three figures may be interpreted as giving the various stabilities 
and instabilities of the system, just as was done in the first case. 


The solution of the problem, which has been given and discussed, 
above, gives merely the sequence of events, and does not show the rate 
at which the changes in the system, take place. It will now be shown 
how the time may be found as a function of x. 

Consider the equation 


dx 

dt 


sin 



/is here the angle of lag of the sidereal semi-diurnal tide of speed 2?i ; 
then by the theory of the tides of a viscous spheroid, tan 2/=2??/p, 
where p is a certain function of the radius of the planet and its density, 
and which varies inversely as the coefficient of viscosity of the 
spheroid.^ 

Since by hypothesis the viscosity is small, / is a small angle, so that 
sin 4f may be taken as equal to 2 tan 2 f. Thus, sin 4//-a is a constant, 
depending on the dimensions, density, and viscosity of the planet. 

It has already been shown that t 3 varies as a?"* 13 , and g is a constant, 
which depends only on the density of the planet. Hence, the above 
equation may be written 

x lz —=K(n-Q), 
dt ’ 


where K is a certain constant, which it is immaterial at present to 
evaluate precisely. 

Since n~h—x and Q=aT 3 , we have 


K dt=z 


— x 15 dx 
x 4 — hx s - hi’ 


or — 1 —-——— -f a const. 

)x 4 -hx 8 + l 

The determination of this integral presents no difficulty, but the 
analytical expression for the result is very long, and it does not at 
present seem worth while to give the result. The actual scale of time 
in years will depend on the value of K, and this is a subject of no 
interest at present. 

It will, how r ever, be possible to give an idea of the rate of change 
of the system without actually performing the integration. This may 


* cc On the Bodily Tides of Viscous and semi-elastic Spheroids,” &c. “ 

Trans,,” Part I, 1879, p. 13, § 5. 


Phil. 
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be done by drawing a curve in which the ordinates are proportional 
to dt/dx, and the abscissae are x. The equation to this curve is then 

■gdt _ —a? 15 

tlx x 4, ~h'P ~{~ 1 

The maximum and minimum values (if any) of dt/dx are given by the 
real roots of the equation 

ll®*~12fo 3 + 15==0. 

One of such roots will be found to be intermediate between a and b, 
and the other greater than a. 


Fro. 8, 



Diagram illustrating the Rate of Change of the System. 


Fig. 8 shows the nature of the curve when drawn with the free 
hand. It was not found possible to draw this figure to scale, because 
when h—2'6 it was found that the minimum M was equal to *85, and 
could not be made distinguishable from a point on the asymptote A, 
whilst the minimum m was equal to about 900,000, and could not be 
made distinguishable from a point on the asymptote C. 

The area intercepted between this curve, the axis of x , and any pair 
of ordinates corresponding to two values of x, will be proportional to 
the time required to pass from the one configuration to the other. 

Where dt/dx is negative, that is to say, when the satellite is falling 
into the planet, the areas fall below the axis of x. This is clearly 
necessary in order to have geometrical continuity in the curve. 

The figure shows that the rate of alteration in the system becomes 
VOL. sxx. u 
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very slow when the satellite is far from the planet; this must indeed 
obviously be the case, because the tidal effects vary as the inverse sixth 
power of the satellite’s mean distance. 


V. “ On the Modifications of the Spectrum of Potassium which 
are Effected by the Presence of Phosphoric Acid, and on 
the Inorganic Bases and Salts which are found in combina¬ 
tion with Educts of the Brain.” By J. L. W. Thudichum, 
M.D., F.R.C.P.L. Communicated by John Simon, C.B.. 
F.R.S., &c. Received March 10, 1880. 

Among the results of a large investigation on which I have for 
many years been engaged in regard of the chemistry of the brain, I 
had been led to conclude that the so-called “protagon” of Oscar 
Liebreich is not a definite chemical body, but is a variable mixture of 
several bodies. This conclusion of mine (which agrees with opinions 
expressed on the same subject by Streeker, Diaconow, and Hoppe- 
Seyler) was published by me in 1874,* and endeavours to controvert 
it have since then been made, on several occasions, by Dr. Arthur 
Gamgee.-f Last summer, he brought before the Royal Society J his 
contentions for the chemical individuality of “protagon”; and it 
fortunately was in my power shortly afterwards to publish evidence, 
which, I believe, those who will take the trouble to follow it will 
find quite unanswerable, that Dr. Gamgee’s contentions were mis¬ 
taken^ Part of my evidence to that effect consisted in showing by 
quantitative analyses that Dr. Gamgee’s so-called “protagon” con¬ 
tains 0*7 per cent, of potassium; secondly, that in connexion with 
trifling differences in the extraction process, the proportion of potas¬ 
sium in different specimens of “ protagon ” can be made to range 
from a trace to 1*6 per cent.; thirdly, that with the variable quantities 
of potassium the quantities of phosphorus and other ingredients will 
also vary. 

In the last published number, TSTo. 200, p. Ill, of the “ Proceedings 
of the Royal Society,” I find that Dr. Gamgee has recently brought 
the question again under notice of the Society, and that, in doing so, 
he especially rests his case upon the following statement made by his 
colleague, Professor Roscoe, on the subject of some examinations, 

# “ Reports of the Medical Officer of the Privy Council and Local Gfovernment 
Board.” New Series. No. III. 

+ £ * Zeitschrift fur Physiol. Chemie,” vol. iii, p. 260 $ a Ber. Deutsch. Chem. 
Gres.,” 1879, &c. 

X “Proc. Roy. Soc.,” vol. xxix, p. 151. 

§ “ Annals of Chemical Medicine.” Edited by J. L. W. Thudichum. Vol. i, 
p. 254. 
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Diagram for Inclination of Satellite’s Orbit.—First ease. 
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Equator.—Second case. 
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Diuirnmi for Inclination of Satellite's Orbit.—Second case. 
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Diuirnmi for Inclination of Satellite's Orbit.—Second case. 
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Diagram illustrating the Rate of Change of the System. 

















